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Abstract 

We examine quark flavour mixing matrices for three and four generations using the recursive 
parametrization of U(n) and SU(n) matrices developed by some of us in Refs. [2] and [3]. After a 
brief summary of the recursive parametrization, we obtain expressions for the independent rephas- 
ing invariants and also the constraints on them that arise from the requirement of mod symmetry 
of the flavour mixing matrix. 
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I. INTRODUCTION 



Study of flavor mixing in weak interactions provides a low energy window for new physics. 
Currently, experiments are underway at Belle and BaBar to check the "unitarity triangle" for 
the 3x3 flavor mixing matrix, as accurately as possible. If there is a significant deviation 
then it would be a signal for the existence of more than three generations. Furthermore, 
the 3x3 CKM mixing matrix contains only one CP-violating phase thus implying that 
CP-violations in different processes are related. Again, the violation of any one of these 
relations would be a signal for more generations. Consequently, in this paper we study some 
general properties of a 4 x 4 flavor mixing matrix. Such a matrix in general has six angles 
and three phases. However, a moduli symmetric 4x4 unitary matrix has fewer parameters. 
We study such a matrix in detail and present parametrizations which would be useful for 
confrontation with experiments in the future. 

In Sec. II, rephasing invariants for a n x n unitary matrices are defined. In addition, 
relations between plaquettes for the particular cases n = 3 and 4 are given. In Sec. Ill, 
recursive parametrization for the n x n case is given together with that for n — 2, 3, and 4. 
Rephasing invariants in the recursive parametrization are presented in Sec. IV. A moduli 
symmetric unitary matrix has fewer parameters and the results for n = 3 and 4 are given 
in Sec. V. In Sec. VI the standard PDG parametrization [1] is obtained using the recursive 
approach. The conclusions are presented in Sec. VII. 

II. REPHASING INVARIANTS OF U(n) MATRICES 

It is known that for a given U{n) matrix (n x n unitary matrix) V under rephasing, i.e., 
under multiplication by independent diagonal U{n) matrices from the left and the right, 

V^V' = D(6')VD(6), (1) 

with 

D(ff) = diag{e*\ . . . , e i6n ) and D{6') = diaglV* 1 ', . . . , e<}, (2) 

the basic quantities that remain invariant are (a) the n 2 moduli \Vij\ = 1, . . . ,n) and 
(b) the n 2 (n — l) 2 /4 Bargmann invariants or "plaquettes", A^i = VikV* k VjiV£ (i < j, 
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k < I; k, I — 1, . . . , n). Unitarity of V gives rise to relations between them and one finds 
that there are n{n — l)/2 independent invariants of type (a), i.e., of modulus-type. As to 
those of type (b), the phase-type invariants, using unitarity one finds that the n 2 (n — l) 2 /4 
invariants can all be algebraically expressed in terms of (n — l) 2 elementary plaquettes 
= Aj j +lj fc 5 fc +1 . The {n — l) 2 basic plaquettes, in turn, are related to each other and can 
be expressed in terms of (n — l)(n — 2)/2 algebraically independent primitives which may be 
taken to be A ik (i < k < n — 1). The total number of algebraically independent rephasing 
invariants is thus n(n — l)/2 + (n — l)(n — 2)/2, i.e., (n — l) 2 . For n = 3 and n = 4, the 
relations are explicitly given below. 

A. Relations between plaquettes for n = 3 

Here there is only one primitive, viz., A i2 . Orthogonality of the rows of V gives the 
relations, 

A 4l = -|\/, 2 | 2 |\/ m2 | 2 -A i2 . (3) 
Likewise the orthogonality of the columns gives, 

A lt = -|\/ 2l | 2 |\/ 2m | 2 -A 2j . (4) 

These are four inhomogeneous equations for four quantities An, A 21 , A 22 , and A 12 . One 
of them is derivable from the other three leaving us with three equations which allow us to 
solve for An, A 21 , and A 22 in terms of A 12 : 



An = 


— ^12 2 ^22 2 — A^ 2 , 




A 22 = 


— V>2 2 2 V 2 3 2 — A^ 2 , 


(5) 


A 21 = 


|^ 22 | 2 (|\/ 23 | 2 -|\/ 32 | 2 ) + A 12 . 




Any other plaquette, e.g., A 1213 , 


can be expressed as A n A 12 / \Vi 2 2 V 22 2 and, using 


the 



relations above, as — | Vi 3 | 2 | V 2 3| 2 — A 12 . 

As is well known, these relations have the consequence that the imaginary parts of all the 
plaquettes are the same, up to a sign. Furthermore, if even one Vij, say Vn, vanishes, then 
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all the plaquettes become real. It is also evident that, imposing mod symmetry on V, i.e., 
requiring \Vij\ = \Vji\, while reducing the number of independent modulus type invariants 
from three to two, has no effect on the number of independent phase type invariants. 



B. Relations between plaquettes for n = 4 

In this case row and column orthogonality of V respectively yield: 

\V i3 \ 2 \V i+13 \ 2 (^n + A* 2 ) = -A* 2 (A, 2 + A* 3 ) (i = 1, 2, 3), (6) 

\V 3l \ 2 \V 3t+1 \ 2 (A ll + A* 2i ) = -A* 2l (A 2t + A* 3l ) (i = 1,2,3), (7) 
which may alternatively be written as 

|^ 2 | 2 |^ +12 | 2 (A, 3 + A* 2 ) = -A* 2 (A i2 + A*J (i = 1,2,3), (8) 

\V 2i \ 2 \V 2i+1 \ 2 (A 3i + A* 2i ) = -A* 2l (A 2t + A* lt ) (i = 1,2,3). (9) 
Choosing Eqs. (6) and (9) we have: 



\v 13 


2 


^23 


2 (A n 


+ a; 2 ) = 


-a; 2 (a 12 


+ Ay, 


\v 23 


2 


^33 


2 (A 21 


+ Ay = 


— A 22 (A 22 


+ Ay, 


\v 33 


2 


v 43 


2 (A 31 


+ Ay = 


-A* 2 (A 32 


+ Ay, 


\v 21 


2 


v 22 


2 (A 31 


+ Ay = 


-A^(A 21 


+ An), 


\v 22 


2 


v 23 


2 (A 32 


+ Ay = 


— A 22 (A 22 


+ Ay, 


\v 23 


2 


v 24 


2 (A 33 


+ Ay = 


-a; 3 (a 23 


+ Ay. 



These six equations for the nine plaquettes allow us to solve all of them in terms of the 
primitives which we choose to be Ai 2 , Ai 3 , and A 23 . The relevant equations are: 

An = -m,pny ' - a; 2 (l + i^jfr^) (ID 
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Aa, = HV„lW-^(l + ? ^ ? ) (12) 

+ dfep) ~ + ^w) A - = ™ a <™ a - '^ |2) (13) 

( 1+ ^) A -( 1+ ^fc) A3 - 

i^i 2 ^ 2 (i + T ^ gF ) - i^n^i 2 (i + d4) 

A 31 = -|V 31 | 2 |V 32 | 2 - A* 21 (l + p^r-p ) (15) 
A 22 = 21 + 1 32 ' ■ (16) 

III. RECURSIVE PARAMETRIZATION OF C7(n) (<Sl7(n)) MATRICES 

Let £/(n) denote the group of unitary matrices acting on all n dimensions. For to = 
1,2,. — 1, we will denote by U (to) the unitary group acting on the first to dimensions, 
leaving the dimensions to + 1, . . . , n, unaffected. Then we have the canonical subgroup chain 

17(1) C 17(2) C ••• C U(n- 1) C U(n). (17) 

General matrices of U(n), U(n — 1), . . . will be written as A n ,A n -i, . . ., respectively. In a 
matrix A m G [/(to) the last rows and columns are trivial, with ones along the diagonals 
and zeros elsewhere (when no confusion is likely to arise, A m will also denote an unbordered 
to x to unitary matrix). 

It was shown in [2] that any matrix A n G U{n) can be expressed uniquely as an n-fold 
product 

An = MOAn-^An-iiZ) • • • Mi)MP)M<*)Mx) , (is) 
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where A n (C) is a special U(n) element determined by an n-component complex unit vec- 
tor £, A n -i(ji) is a special U(n — 1) element determined by an n — 1-component com- 
plex unit vector rj, and so on down to .4.2 (cO that is a special U{2) element determined 
by a two-component complex unit vector a, and A±(x) is a phase factor belonging to 
U(l). The complex unit vectors {£,77,...}, appear as the last columns of the (unbor- 
dered) matrices {A n (C), An-iivi), . . .} and can be identified with the labels of the cosets 
{U(n)/U(n — 1), U(n — 1)/U(n — 2), . . .}. Remembering that {£, 77, . . .} are complex unit 
vectors of dimensions {n,n — 1, . . .}, it is easily seen that the number of real independent 
parameters add up to n 2 as they should. 

The same considerations as above apply to SU(n) matrices as well. Denoting by A n (£) 
the corresponding matrices in SU(n), any A n e SU(n) can be decomposed as 

A n = A n (C)A n _ 1 (r 7 )A n _ 2 (e) • • • M~/)MP)M<*) (19) 

The above construction fixes only the last column of the unitary matrix A n (C) as C 
and one has a great deal of freedom in arranging the remaining n — 1 columns leading to 
many explicit forms for these matrices. In this work we consider two explicit forms which 
correspond to those discussed in [2] and [3] respectively. 

The explicit expressions for the nonzero matrix elements of A n (C) considered in [2] are 



A n (C) = (a jk (C)) e U(n); j, k = 1, 2, . . . , n. 

MC) = 0; J = l,2,...,n. (20) 
a ii-i(C) = Pj-i/pj', J = 2, 3, . . . , n; 

P 3 = \/l ~ IC J + l| 2 " IC, +2 | 2 ICnl 2 = ^Kll 2 + • • ' + I0I 2 . 

%fe(0 = -QCk+i/PkPk+i, j < k < n - 1. 
Thus, for instance, for n = 2, 3, 4 we have: 

^l 2 (a) =| I , /ii = (21) 
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/ -p*M ai a 2 -FM<t 2 fa ^ 



AM 



0-1/0-2 





-PZfafa fa 

0"2 fa 



= \fa\, °i = ViiWTW- (22) 



Mi) 



( -lll\lp\p2 ~lsll/p2p3 ~llll/p3 ll\ 

Pi/ Pi -itli/piPz -llli/pz 72 
P2/P3 -1II3/P3 73 

p 3 74 / 

P3 = VlTll + |72| + |73| 2 - 



V 

pi = M, P2 = Vliil 2 + 72I 2 , 



(23) 



The determinant of the matrices vA n (£) turns out to be (— l) n 1 Ci/|Ci| an d hence the 
corresponding SU (n) matrices can be obtained by multiplying, for instance, the first column 
by (-I) n-1 G7IG|- Thus for n = 2,3,4 we have 



A 2 (a) = 




(24) 



/ _p*/ a2 -fa 3 fa/a 2 fa ^ 



AM 



\ 



fal/°2 





(25) 



/ 



^4(7) = 



(26) 



-fa 3 fa/*2 fa 

o- 2 fa 

I l 2 /p2 -ltll/p2pZ ~llll/p3 ll\ 

-ill Pi -I3I2/P2P3 -1II2/P3 72 
Pi/p 3 -7473/P3 73 

p 3 74 / 

This parametrization assumes that Ci is nonzero. As a result, in the extreme case when 
C = (1, 0, . . . , 0), the matrix A n (C) does not reduce to the identity matrix. A parametrization 
where this does happen and which corresponds to that given in [3] is given below: 



V 
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A n (C) = (a jk (C))eU(n); j,k — 1,2, ... ,n. 

a jn (C) = 0; j = l,2,...,n. 

a jj(C) = Pj/pj-i, j = 1,2,..., n- 1; p = 1; 

p 3 = 7 1 -Kil 2 -IC2l 2 I0l 2 = V / IO + il 2 + 

a ife(C) = -QCk/iPk-iPk); 3>k, k = 1,2,..., n- 1. 



(27) 



Note that we are using the same symbols as in the parametrization earlier though with 
different meanings. For n — 2, 3, 4, we explicitly have: 



As(/3) 



«4 2 (a) 



Pi «i 
-ala 2 /pi ot2 

\ 



Pi = \ot 2 \. 



<?1 ft 

-PWl/Vl C2/C1 #2 



0"! 



(28) 



0"2 



(29) 



A(7) = 



/ Pi 

-7i72/>i P2/P1 
-ills/ Pi -72 7s/piP2 P3/P2 



7i 
72 

73 



\ 



\-7i74/Pi -7274/P1P2 -7I74/P2P3 74/ 



Pi = VI72I 2 + |7s| 2 + |74| 2 , P2 = VI73I 2 + |74| 2 , P3 = M- 



(30) 



The determinant of the matrices A n (C) is Cn/|Cn|- We can convert the above matrices to 
SU(n) matrices by multiplying, say the (n — l)-th column by C/ICn|- Thus, for n = 2,3,4, 
we have: 



A 2 (a) 



a 2 di 



MP) = 



01 (3 1 

-fifo/<Tl &/<Tl (32 

{ -m/n -ft/a & j 



(31) 



(32) 
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Ml) 



(33) 



(pi 7i \ 

-I1I2/ Pi P2/P1 72 

-ills/ Pi -I2I3/P1P2 ll/p2 73 

\-7i*74/Pi -7274/P1P2 -73/P2 74/ 
Given a matrix *4„ e U(n), we can determine the parameters, the complex unit vectors, 
{£, 77, . . .} in a recursive fashion through the following steps. 



Write A n = {(ijk) £ U (n) as 



»4-7i A n {C^)18 n — 1, 



(34) 



where £ is the last column of „4„ 



~~ a j«- 



(35) 



With A n (C) thus determined, we have 



(36) 



The matrix elements (fry), n — 1 > i , j > 1 , of £>„_i in the first form [2] are given by 



n 1 " 
% = ^2 a li(C) a kj = — ^2 a ln( a kj a i+ln - akndi+lj) 



k=l 



PiPi+l 



(37) 



k=i+l 



and in the second form [3] by 



n 1 n 

bij = ^ ] a ki(C) a kj = ^ ] a kn( a ij a kn ~ a in a kj) ■ 



k=l 



PiPi-1 



(38) 



k=i+l 



Write B n _i as 



# n _i = A n -i{j])C n - 2 



(39) 



with 
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Vj = bjn-i- (40) 
• Repeat the same procedure as above with C n - 2 . 

The same procedure applies to the decomposition of an SU (n) matrix. Thus, for instance, 
using the second form [3], a matrix V G SU(3) can be decomposed as 

V = A 3 ((3)A 2 (cx), (41) 

where 

Pi = V l3 , (3 2 = V 23 , (3 3 = V 33 , (42) 

and 



[V 2 * 3 (V 12 V 23 - V 13 V 22 ) + V 3 * 3 (V 12 V 33 - V 13 V 32 )] (V 33 V 22 - v 32 v, 



L ' 1'6 V ' J-^ ' ^ ' ±o r zzy i * * * oo ' lo ' oz / j v ' oo • zz • oz • 23 ) / . \ 

v^W+F^F VI^3| 2 + |^ 33 | 2 

IV. REPHASING INVARIANTS IN THE RECURSIVE PARAMETRIZATION 

Having shown how to parametrize a given U{n) (SU(n)) matrix in terms of a sequence 
of complex unit vectors {£, rj, . . .} of dimensions {n, n — 1, . . .}, we now examine how these 
parameters transform under rephasing with the purpose of constructing rephasing invariants 
out of them. For simplicity and without any loss of generality we will assume that the 
given matrix V belongs to SU(n) and will consider the cases n = 3,4 and discuss the 
transformation properties of the parameters in both the forms [2, 3] given above. In the first 
form [2] , any SU (3) can be written as 

V = A 3 ((3)A 2 (<x) 



( -/3%a* 2 /cr 2 + /^/3ic^/<7 2 ~P 2 oii - /3 3 /3ia 2 a 2 ft ^ 



(44) 



f3{a* 2 /a 2 + Pl(5 2 a\/a 2 (3{ai/a 2 - /3 3 /3 2 a 2 /a 2 f3 2 
y -cr 2 al a 2 a 2 ft J 

Under rephasing by independent diagonal SU(3) matrices D(6) and D(9') where D(9) = 
diag(e^ ei+d2 \ e i( - 9l+d2 \ e~ 2id2 ) and D(9') is similarly defined, we have 
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V ->• V = D(9')VD(9) 

= A 3 (f3')A 2 (<x') (45) 

From the locations of a±, a 2 , Pi, (3 3 , and (5 3 in Eq. (44) one can easily deduce the transfor- 
mation properties of j3 and a: 

a[ = a^- 6 ^, a' 2 = a 2 e^- e ^\ (46) 



f3[ = f3 ie ^ +e '^\ (3' 2 = I3 2 e l{ - e '^- 2e >\ (3' 3 = f3 3 e 1 ^- 26 ^. (47) 

From these transformation properties it is evident that (aia^PlP^Ps) an< ^ hence 
arg (aia 2 (3l(3 2 (3 3 ) is invariant under rephasing. 
For n — 4, parametrizing D{9) as 

D{9) = dmg{e i ^ +e ^\e ti ^ +9 ^\ e ~ 2ie * +ie \ e" 3 ^ 3 ) (48) 
and similarly for D(9'), one finds that 



V = A 4 (^)A 3 (f3)A 2 ( a ) -> V = D{9')VD{9) 

= A 4 ( 7 >)A 3 ((3>)A 2 (cx') 

= D(9')A 4 {~f)di&g{e W3 , e ie \e ie \ e -3 ** 3 ) 

A3(^)^ 2 (a)diag(e iei+lf?2 , e^ 1 ^ 2 , e~ 2ie \ 1) (49) 

The expressions for y can easily be read off: 



sy[ = 7ie *(*i+*a+«3-3«3) ) y = 7ie i(-ei+^+^3-3»3) > 

y 3 = 7ie *(-^3-3«s) > y = 7ie -si(fl'3+3«3) - (50) 
A little algebra shows that 



D(^)^4(7)diag(e if?3 , e ifl3 , e ie \ e" 3 ** 3 ) = 

A 4 (y)diag(e^ 2 ^ +29 ^ 2e 3) ? e i(-2^+flS+«3) > e <(-3«S+fl3) > 1)? (51) 
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so that the rest reduces to an SU (3) problem in a 3 x 3 matrix form 



MP)M<*') = diag(e i ^ +2e 3- 2e 3) ) ^(-2^+^+e3) )e i(-3^ + e 3 ) ) 1} 
A 3 (f3)A 2 (a)dia,g(e Wl+l0 \ e~ m+ie \ e" 2 ^ 2 ). 



(52) 



We see that for the SU (4) problem to accompany Eqs. (50) we have, 



a 1 = a±e 



,1(30^-61-02-03) 



a 2 = a 2 e 



,i(-30^-0i+0 2 +0 3 ) 



(53) 



ft = ^(202+203-202-203) ; ^ = ^ 2 gi(-20 2 +0 3 -202+03) ^ ^ = ^ 3 gi(-30 3 -202+03) _ (54) 



With the transformation properties of 7, (3, and a at hand, we can now systematically 
construct rephasing invariant quantities out of them as shown in Ref. [2]. The three indepen- 
dent invariants turn out to be (a^a^*/^/^), (P2P3I3I4), and (Pi^lil^ls)- The arguments 
of these quantities furnish the three independent phase type invariants for the 577(3) prob- 
lem. Notice that the first of these is the rephasing invariant for the S77(3) problem and this 
is indeed a rather desirable feature of the recursive parametrization outlined here as one 
goes from n to n + 1 one retains the parameters at the n th level. 

In the second form [3], for n — 3, the analogues of Eqs. (44), (46), and (47) are 



V = A 3 ((3)A 2 (a) 




(55) 



a 1 = a±e 



,i(0£ +0 2 -01 +0 2 ) 



a 2 = a 2 e 



i(0^+0 2 +0i+0 2 ) 



(56) 



-2i(0 2 +0 2 ) 



(57) 



and the rephasing invariant is (aia 2 (3l(3 2 (3 3 ). 

For n = 4 the corresponding equations to (50), (53), and (54) are 
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sy[ = 7ie ^+ e 2+ 9 3- 393 ), 



7 3 



7ie 



i(-20 2 +0 3 -30 3 ) 



y 2 = 7l e i (-«i+^-3«3) > 
y = 7l e- 3 ^a +3f? 3). 



(58) 



aie 



i(6>i+^+^-6»i+6l2+6»3) 



a 2 = a2e -<(ei+^+^3+»i+« 2 +«3). 



(59) 



(60) 



The three independent invariants turn out to be (aia^/?^!/^), (PiP^li^), an d 

(^2^72*7374) ■ 

V. CONSTRAINTS DUE TO MOD SYMMETRY 

In this section we examine the constraints on the parameters that arise from demanding 
that the given SU(n) matrix be mod symmetric, i.e., \Vij\ = \Vji\. For convenience we shall 
use the first form [2] for this discussion. 

For n — 3, mod symmetry requires that \a 2 \ = |ft|/c2- The number of independent angle 
type invariants comes down from three to two leaving the phase type invariant unchanged. 

For n = 4, with V = At(7)A 3 (/3)A 2 (a:), after some algebra one finds, 



(61) 
(62) 



IK4l = 


v 41 \ 


=>• |at 2 | = 72I/P2, 


\Vu\ = 


v i3 \ 


=> |ft| = I73IM, 


\v 23 \ = 


v 32 \ 


=>• cos((<Ji + 5 2 + $0/2) x 




rift 


-cos((5 1 - 5 2 - 5 3 )/2) + ^ 








denote 


the 


three independent invariant 



= 0. 



(63) 



arg(/3 2 /3|7^7 4 ), and arg(/3!/3 2 * 7^7^73), respectively. The equalities \V 24 \ = \V 42 \, \V 12 \ = \V 21 \, 
and IV13I = I V31 1 give no new conditions. It can be seen from the above equations that one 
can obtain mod symmetry by requiring 



\a 2 \ = I72I/P2, |ft| = I73I/P3, S 1 + 5 2 + 5 3 = 7T, 



(64) 
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and in this situation the mod symmetric matrix mixing is parametrized by four angles and 
two phases. 

A simpler moduli symmetric parametrization can be obtained if some of the eigenvalues 
Ei (i = 1,2, .. . ,n), of the n x n unitary matrix V are equal. For the case when n — 1 
eigenvalues are equal, viz., E 2 — E 3 — • ■ • — E n , V can be expressed in terms of n — 1 real 
parameters and only one phase [4]. 

VI. COMPARISON WITH THE "STANDARD" (PDG) PARAMETRIZATION 

For the case of three generations, the standard or PDG [1] parametrization of the mixing 
matrix is obtained by putting 

oil = ci 2 , ol 2 = S12, (3i = s 13 e~ lSl3 , (3 2 = S23C13, As = C23C13, (65) 

(cy = cosOij and = sin^) in V = A 3 ((3)A 2 (ol) with A 3 ((3) and A 2 (a) given by Eqs. (32) 
and (31), respectively. 

The extension of the mixing matrix to four generations is given by V = A^'y) A 3 (f3) A 2 (ct) , 
where ^4(7) is given by Eq. (33) with (3 and a as before and 

7i = si4e~ h \ 72 = c 14 s 2 ^ l5 ' M , 73 = S34C24C14, 74 = C34C24C14, (66) 

which conveniently reduces to the case of three generations when 9 U , 9 24: , and #34 are all set 
equal to zero. 

We note here that the parametrization given above is closely related to the Harari-Leurer 
parametrization [5] where the mixing matrix is expressed as an ordered product of essentially 
2x2 "rotation" matrices. Our parametrization results when one suitably combines the factors 
appearing in that form. For instance, in the 4x4 case, the Harari-Leurer form for the mixing 
matrix has the structure ^234^24^14^23^13^12 and reduces to our form by the identifications 
^4(7) = ^34^24^14, A 3 (f3) = fi 2 3^i3, and A 2 (ac) = Qi 2 , provided we choose 74, (3 3 , and a 2 
to be real. 
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VII. CONCLUSIONS 



In this work we have examined in detail the question of parametrizing quark flavor mixing 
matrices for three and four flavors within the framework of the recursive parametrization 
developed in Refs. [2] and [3]. In particular we have shown, given the matrix, how to 
determine the corresponding parameters. We have also studied in detail aspects of rephasing 
invariants in this parametrization scheme and have derived conditions for the mixing matrix 
to be moduli symmetric. 
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